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Abstract 



We calculate the effects of vertex corrections, of non-constant density of states and 
of a (self-consistently determined) phonon self-energy for the Holstein model on a 
3D cubic lattice. We replace vertex corrections with a Coulomb pseudopotential, 
adjusted to give the same Tc, and repeat the calculations, to see which effects 
are a distinct feature of vertex corrections. This allows us to determine directly 
observable effects of vertex corrections on a variety of thermodynamic properties of 
superconductors. To this end, we employ conserving approximations (in the local 
approximation) to calculate the superconducting critical temperatures, isotope co- 
efficients, superconducting gaps, free-energy differences and thermodynamic critical 
fields for a range of parameters. We find that the dressed value of A is significantly 
larger than the bare value. While vertex corrections can cause significant changes in 
all the above quantities (even when the bare electron-phonon coupling is small), the 
changes can usually be well-modeled by an appropriate Coulomb pseudopotential. 
The isotope coefficient proves to be the quantity that most clearly shows effects of 
vertex corrections that can not be mimicked by a /x^. 
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I. INTRODUCTION 



The theory of conventional, low-temperature superconductors has been well understood for 
decades, within BCS theoryi, and many material properties of superconductors have been accu- 
rately described within the more appropriate Migdal-Eliashberg formalismSa, which includes the 
retardation effects of the electron-phonon interaction in a realistic manner. The success of the 
formalism arises from Migdal's theoremi, which says that vertex corrections can be neglected when 
the ratio of the phonon energy scale to the electron energy scale is small (such as the value of 10~^ 
typical of conventional low-temperature superconductors). The physical reason being that the ion 
movement is typically too slow to respond to anything but the mean-field potential produced by 
the fast-moving electrons. 

In recent years there has been investigation of superconducting materials!^ such as 
Bai_2^K^Bi03, K^Cqq and the Al5s, whose phonon energy scale is a larger fraction of their elec- 
tron energy scale. For such materials, the standard Migdal-Eliashberg theory may no longer be 
valid, either because second-order orocesses in the electron-phonon coupling (so-called vertex cor- 
rections) can not be neglectedEJJTJ, or because structure in the electronic density of states and 
finite-bandwidth effects become significantEl. Our aim is to identify those experiments which can 
clearly indicate where vertex corrections, or structure in the electronic density of states, are observ- 
able in real materials. In particular, we wish to uncover those effects which can not be mimicked 
by a Coulomb pseudopotential, /x^, as this would render them unobservable in practice, because a 
IIq is typically fitted to the experimental data. 

The most compelling evidence for the effect of vertex corrections would come from a tunneling 
conductance measurement that provided data for the tunneling conductance out to an energy at least 
twice that of the maximal phonon frequency of the bulk material. Then a tunneling inversion could 
be performed using only the experimental data out to an energy of the maximal bulk phonon energy, 
and the results for Migdal-Eliashberg theory would be compared to the vertex-corrected theory for 
the experimental data that was measured in the multiphonon region, at voltages above the maximal 
bulk phonon energy of the material. Unfortunately, such an analysis has only been performed for 
leaded and in that case, the experimental data was not accurate enough in the multiphonon region 
to be able to see if effects of vertex corrections were observable. Similarly, high energy data from 
optical conductivity experiments would indicate whether or not vertex corrections are important, 
but there again, the accuracy of the data might preclude seeing effects of vertex corrections. 

In order to determine what effects are unmistakenly due to vertex corrections, we fit a Coulomb 
pseudopotential to a Migdal-Eliashberg theory so that the superconducting transition temper- 
ature Tc is the same for the vertex-corrected theory (with no /i^) and the conventional Migdal- 
Eliashberg theory (with a /i^ but no vertex corrections). In the case of dressed phonons, we adjust 
both the electron-phonon interaction energy and /i^, in order to fit the same values of A, and as 
the results with vertex corrections. This procedure is exactly the procedure carried out in analyzing 
experimental data within the conventional theory. We carry out this "experimental analysis" on a 
model system where we know that vertex corrections have a large effect, in order to see the extent to 
which the conventional analysis masks their observation. Hence we look for differences in thermo- 
dynamic quantities between these results and those with vertex corrections. Any quantities which 
are significantly different when calculated with a fitted fi^ instead of vertex corrections pin-point 
the experiments that can give the best indication that a material needs to have vertex corrections 
included in its description. 

Since this is the initial attempt at solving such a problem, we study a simple model system, which 
has a non-constant electronic density of states (Figure |lD, using both Migdal-Eliashberg formalism. 
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and going beyond it, to include vertex corrections. The vertex corrections are second-order diagrams 
where a pair of phonon hnes cross, as shown in the self-energies of Figure 0(b). In both cases, the 
phonon propagators can be dressed {i.e. with a phonon self-energy included) or bare, as shown 
in Figure |^. We find that dressing the phonons leads to a strong enhancement of the effective 
interaction strength, exemplified by a large renormalization in the value of the electron-phonon 
coupling parameter, A . 

The specific model we study is the Holstein HamiltonianO on a 3D cubic lattice. The Hamil- 
tonian consists of conduction electrons that hop from site to site, coupled to harmonic, localized 
(Einstein) phonons: 



where cfg. and Qo- are fermionic operators which create and destroy, respectively, an electron of 
spin-o" in a single Wannier (tight-binding) state on the lattice site ^, whose total electron occupancy 
is given by rij = cJt^iT + The electron hopping is between nearest neighbors only, such that 

tij = —t if i,j are neighboring sites, with t the overlap integral, and tij = otherwise. The phonons, 
of mass, M, with displacement Xi and momentum, pi, are characterized by their frequency, Q. The 
strength of the bare electron-phonon coupling can be measured by the bipolaron binding energy, 
U, where: 

The chemical potential is /i, and is always calculated self-consistently for a given average number 
of electrons per site, n (0 < n < 2). Particle- hole symmetry occurs at half- filling, where n = 1 and 
H = U. We concentrate our work on superconductivity and ignore any possible charge- density wave 
order that may occur near half-filling. 

We carrvout weak-coupling expansions within the conserving approximations of Baym and 
KadanoffE^iEl. The electron self-energy, S, is given as a functional derivative of the free-energy 
functional, $, with respect to the electron Green's function, G. When dressing the phonons, we 
maintain a conserving approximation by careful choice of the phonon self-energy. In this case, 
there must be a free-energy functional whose partial derivative with respect to the electron Green's 
function yields the electron self-energy, and whose partial derivative with respect to the phonon 
Green's function, D, yields the phonon self-energy. In such cases there always exists a related 
free-energy functional whose full functional derivative with respect to the electron Green's function 
yields the electron self-energy. In all cases, the irreducible vertex function is given by the full 
functional derivative of the electron self-energy with respect to the electron Green's function, i.e. 
T = 5T./5G. 

We employ the local approximation in our calculations, which means we neglect the momentum 
dependence in the self-energy and irreducible vertex functions. The local problem gives an exact 
solution in the infinite-dimensional limitlllll, but in this case it is an approximation, which, as used 
in Migdal-Eliashberg theory, is expected to give good quantitative results, though of course it means 
we can only study s-wave pairing. 

We use the formalism, described in detail in the next section, to calculate the critical tempera- 
ture, Tc and the isotope coefficient, a, by looking at the instability of the normal state. Also, within 
the superconducting state, we calculate the superconducting gap. A, and the thermodynamic crit- 
ical field. He, from the free-energy difference, Fs — F^. We point out in particular those results 
which deviate from Migdal-Eliashberg theory, and arise from vertex corrections or a non-constant 
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density of states. By including a Coulomb pseudopotential, which causes the same change in 
as vertex corrections, we are able to demonstrate what experiments can be used to differentiate 
between vertex corrections and a Coulomb pseudopotential. 

Section II contains the formalism, describing the particular approximations we use, and explain- 
ing how our calculations are carried out. Section III presents our computational results, including 
a number of graphs depicting the various quantities that can be determined experimentally for real 
systems. We make our conclusions in Section IV, which is followed by an Appendix which gives the 
detailed formula withheld from Section II. 



II. FORMALISM 

In our calculations, we use four different types of conserving approximations, which will be 
expounded below. Four different approximations are necessary, in order to reveal how dressing 
the 'bare' phonons affects solutions of the model, as well as to demonstrate the effects of vertex 
corrections. Each approximation includes a specific self-energy, and hence a specific vertex function, 
as shown in Figures In short, these can be described as (a) Migdal-Eliashberg approximation 
with a truncated dressing of phonons; (b) second-order perturbation theory, which contains vertex 
corrections and a truncated dressing of phonons; (c) Migdal-Eliashberg theory with dressed phonons; 
and (d) vertex-corrected theory with dressed phonons. 

The calculations are carried out on the imaginary axisil, with the Green's functions and self- 
energies defined at Matsubara frequencies, iujn = (2?t, -|- l)7riT, where T is the temperature. The 
self-energies are calculated from the derivative of an appropriate free-energy functional, $, as: 
Y^{iujn) = -f ^Q-u^i^ y = sF*^iuj ) ' ^^ilc the irreducible vertex function for superconductivity 

is given by: ^ = sP(Z"\ • Figure ^ shows the Feynman diagrams which correspond to the 
appropriate free-energy functional, $, for each of the four calculations, while we defer the specific 
formulae to the appendix. In the diagrams, taking a functional derivative corresponds to cutting 
and removing the appropriate Green's function line. 

When calculating properties for the superconducting state, we must use the Nambu formalismll, 
where the electron Green's function and self-energy are represented by a 2 x 2 matrix: 

where p''°^(e) is the approximate form for the non-interacting electron density of states on a 3D 
cubic lattice as shown in Figure |l], given by Uhrigil, and 

mcj ) = ( ^^'^"^ '^^'''"^ ] r = ^ ] U) 



Here, the diagonal and off-diagonal Green's functions are defined respectively as: 

G"^'^(z^„) = -^E [dreM^^nr) (T.c,.(r)4(0)) , (5) 

F{iujn) = -j^Y. dTexp{iuJr,T) (T^cj|(r)9|(0)) , (6) 

where T,- denotes time-ordering in r. The definitions of "Eliun) and for each calculation are 

given in the appendix, Eqs. ([A^)-( |AT6|) , and are represented in Figure 0. 
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The phonon propagator, D{iujy), is defined byn3: 

Diiujy) = ——y^ (irexp[it^i,r] (Txj(r)xj(O)) , (7) 
N ^ Jo 

such that the 'bare' propagator, with no self-energy, D^^^icu,^), is equal to ^^^^2 ■ The appropriate 
Matsubara frequencies, u^, in this case, are those which lead to bosonic statistics, such that iuu = 
27riTz/. 

When we use bare phonons, and have no vertex corrections, the electron self-energy is illustrated 
in Figure 0(a). The first term is the Hartree contribution, Un. It can be included by a shift of the 
chemical potential as it only contributes a constant to the diagonal part of the self-energy, Tj{iuJn)- 
Such a shift is implicitly included when fi appears in a Green's function, so that the Hartree term 
is hereafter neglected. It is followed by the Fock term, then a single-phonon dressing term, where 
the phonon line includes a single loop, which is the electron polarizability, 7i^^\iujy): 
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with the factor of two arising from the summation over spin. The third term includes dressing of 
phonons in a truncated manner, and allows us to make comparison with the complete second-order 
approximation, in Figure @(b), where the only difference is the inclusion of the vertex-correction 
term. The specific formula for the extra term in Figure 0(b), which enters our calculations when 
we include vertex corrections using bare phonons, is given in Eqs. ( [A7|) and ([A^). The extra 



contributions to the self-energy coming from vertex corrections have opposite sign to the Fock 
term, near half-filling, where the product of two electronic Green's functions is negative (the Green's 
functions are pure imaginary at half- filling) . Away from half- filling, the Green's functions gain real 
parts, which means that near the band edges the product of two Green's functions can be positive, 
and the vertex-correction terms then add to the Fock term. 

In the calculations with dressed phonons, we have a Dyson's equation for the phonon propagator, 
as depicted in Figure ^ 

D{iuo,) = D^^\iuj,) + D^''\iu,)Yi{iuo,)D{iu,), (9) 

where 'n.{iujjj) is the phonon self-energy. Note that within the conserving approximations, we must 
determine the phonon self-energy by differentiating the free-energy functional: 

m^u) = ^TJ^y (10) 

With no vertex corrections, the phonon self-energy is simply given by the electron polarization, as 
depicted in Figure |^(a). That is, 

Il{iuJ^) = U-K^^\iuJ^). (11) 

In this case, we use the contributions to the electron self-energy that are shown in Figure 0(c), whose 
explicit formula is given in Eq. ([AllI). Note that the term with the single polarization bubble is 



missing, as all orders of such 'necklace' diagrams are included in the Fock term with the dressed 
phonon propagator. 

When we include vertex corrections, as well as dressed phonons, the phonon self-energy gains 
the extra term shown in Figure |5(b). The full expression is given in the appendix [Eq. (|A16|) ]. 
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Note that a fully dressed phonon propagator is included in the phonon self-energy. As shown in 
Figure 0(d), the electron self-energy now has the expected extra term with a crossing of phonon 
lines. It is almost identical to the extra term in Figure |^(b), except of course, now the phonons are 
dressed. Again, the details of the formula can be found in the appendix. 

Our calculations all involve iteration of the Green's functions and self-energies, until a self- 
consistent solution is reached. We begin with the non-interacting Green's functions (set the self- 
energies to zero), and use it to calculate an initial estimate of the self-energies, according to Eqs. ( |A^ 
A16D . The new self-energies are used to calculate updated Green's functions, according to Eqs. (|) 
and (y). The procedure is iterated, so that at each step there is an updated self-energy, which 
includes a fraction of the previous self-energy, the exact fraction variable, dependent upon the 
progress of the iteration. We stop the process when the change in all the self-energies is less than 
one part in 10~^°, which is typically after tens, but sometimes after hundreds of iteration steps. 

Superconductivity occurs below a critical temperature, Tc, where the normal state becomes 
unstable to fluctuations in the pairing potential (the Cooper instability). The instability shows 
itself as a divergence in the pairing susceptibility, Xm,n, which is given by: 

Xm,n = '^rn ^m,n ~~ ^ -^rrt'^r^ (12) 

I 

where F^.n is the irreducible vertex function, to be defined shortly. The bare susceptibility in the 
superconducting channel for momentum, q, is defined as: 

xS^(q) = ^^G(«a;„,k)G(-«cu^,-k + q), (13) 

which becomes in the local approximation (for the zero-momentum pair): 

(o)_Jm[G(«^ 

where 



UJmZ{iuJm) =^m- Im'^{iuJm)\. (15) 

The transition temperature, occurs when the largest eigenvalue of the matrix —TXmTm,i passes 
through unity. We calculate the highest eigenvalue using the power method. 

The irreducible vertex function, Tm,n is given by the functional derivative of the off-diagonal 
self-energy, with respect to the off-diagonal Green's function: 

= 6(f){iuJn)/SF{iuJm). (16) 

As the pairing fluctuations lead to an instability of the normal state, Tn^m must be calculated in 
the normal state [i.e. in the limit F{iujn) i— > 0). 

Figure ^ shows the contributions to the irreducible vertex function for each of the four ap- 
proximations. Each diagram is achieved by removal of one electron Green's function line, from a 
self-energy diagram in Figure 0. The algebraic expressions for each of the four sets of diagrams are 
given in the appendix. 

As we wish to uncover more than the phase diagram given by the different perturbation ap- 
proximations, we move on to describe how we calculate other properties. In order to find the 
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superconducting gap and thermodynamic quantities, calculations are required within the supercon- 
ducting state, but a simple addition to the previous computations in the normal state allows us to 
calculate the isotope coefficient. 

The isotope coefficient, a, describes how the critical temperature changes with the phonon mass, 
M. It is defined as: 

" = -rfh^' 

so that Tc oc M~". The weak-coupling limit of BCS theory, and Migdal-Eliashberg theory with 
no Coulomb repulsion, predict a = 0.5, which corresponds to Tc oc l/\fM. The phonon frequency 
changes with mass, according to 1] oc 1/-\/M, so both the product, MVt^ , and the interaction 
energy, [/, remain constant. Hence we calculate the isotope coefficient simply by changing the 
phonon frequency by 1% (corresponding to a typical mass change of 2% between isotopes) and 
comparing the change in critical temperature. To be precise, we compute a by: 

rp{new) rp{old) ^(old) 

a = 0.5-^ ^j;^ — r — ^, (18) 

rp(old) ^(neui) ^(old) ' ^ ' 

SO the BCS result is achieved if Tc oc 11 

According to standard methodsHIrEa, the energy gap in the superconducting state. A, requires 
a self-consistent calculation within the superconducting state. Note, the order parameter on the 
imaginary axis is related to the off-diagonal self-energy through: 

A{iun) = 4>{iuJn)/Z{iuJn), (19) 

where Z{iujn) is the mass-enhancement parameter, calculated from the electronic self-energy, S(iu;„), 
as given in Eq. ([151) . The gap itself is found from the order parameter on the real axis, at the point 
where Re [A{uj)] = u. We carry out a Fade analytic continuationSill to obtain the order parameter 
on the real axis, and hence the value of the gap. 
The free energy can be found from the formulal22 

F = -2r^ |-ln [-l/det£(zcj„)] + -Trg(icj„)£(ic<;„) j (20) 
„ 1 2 2 J 

+ 1 {In [-1/D{iuj,)] + U{iu,)D{iu,)} + $ + ^(n - 1). 

The free-energy functional, $, is the skeleton-diagram expansion whose differential with respect to 
the electron (or phonon) Green's function gives the electron (or phonon) self-energy [see Figure ^ 
and Eqs. (^), (|A6|), ( [All| ), and ( |A15| )]. We are interested in the free-energy difference between 
normal and superconducting states at fixed electron filling, n, in which case the first, Hartree, term 
in $ is neglected as it is a constant, Un^/2. The final term, fi{n — l), can not be neglected, because 
the chemical potential can differ considerably between the superconducting and normal state when 
one includes the effects of non-constant density of states. 

We calculate the thermodynamic critical field in the superconducting state, H^, from the 
free-energy difference between the superconducting and normal state, according to the formula 
Fs — -Fjv = —fJ-oH^- The thermodynamic field varies with temperature in an almost quadratic man- 
ner, so that calculation of the deviation function, which is defined as the difference between Hc{T) 
and the quadratic form, Hc{0) [1 — (T/Tc)^] gives a sensitive test of changes in thermodynamic 
quantities. 
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In the calculations which include a Coulomb repulsion term Uc, we make the standard simplifi- 
cationli5lE3i23 of only including its effects on the off-diagonal self-energy. This simplification is valid, 
as the normal-state Green's functions in reality include the Coulomb repulsion effects, and these 
change by very little for the diagonal part of the Green's function when superconducting order is 
present. As our model does not include the effects of Uc on the diagonal Green's functions, it is 
not strictly the solution of a simple Hamiltonian with a Uc term includedil. However, the simplifi- 
cation does allow us to compare the effects of including a Coulomb repulsion versus adding vertex 
corrections, on the superconducting properties and transition temperature from a similar normal 
state and it is precisely the method employed in analyzing experimental data on real materials. 
With these comments understood, the only changes to the calculations that are necessary with the 
inclusion of a Coulomb term, are that both the off-diagonal self-energy, (f){iuJn), and the irreducible 
vertex function, Tn,m gain an extra term: 

(l){iujn) ^ (t){iuj„) + UcT^F{iLJm), (21) 

m 

^n,m ^ Tn.m + Uc- (22) 

In the computational calculations, the Matsubara frequency sum is cut off after a constant number, 
Nc, of terms. This leads to a renormalization which reduces the Coulomb term@@ to a pseudopo- 
tential, Uc, given by: 

/ [ 7V,+1 '^rn J 

where the diagonal Green's function, G{iuJm), is at high frequency, where the self-energies may 
be neglected, but must include the self-consistent chemical potential, /x. In our calculations, with 
the frequency cut-off on the scale of the bandwidth such that 256 or 512 Matsubara frequencies 
are used, there is a very small (typically 1% or 2%) reduction in Uc- This is different from the 
conventional approach in real materials because here our energy cutoff is governed by the electronic 
bandwidth, not some multiple of the maximum phonon frequency. To make contact with the 
standard formalism, we define a dimensionless pseudopotential, /i^ = p^^\n) ■ Uc, where p'^'^\jj) is 
the non-interacting electronic density of states at the chemical potential. In calculations at different 
temperatures, Uc is kept fixed, so that /i^ varies to a small extent. 

Finally, we wish to make clear how A, the measure of the electron-phonon coupling strength, is 
defined in our work. A precise definition is required, because different methods of calculating A lead 
to different results away from the weak-coupling limit, especially when the phonons are dressed. 
Here, A is given by: 

A = pW(/i).[/.Z}(0) (24) 

where -D(O) is the zero-frequency component of the dress ec? phonon propagator, which can be sig- 
nificantly different from that of the 'bare' propagator, D^^\Qi). This value of A is usually different 
from the electronic mass enhancement parameter, i.e. A 7^ ^(0) — 1, as the two are only equal in the 
weak-coupling limit and with i-^ 0. The definition of A is identical to that commonly calculated 
from the electron-phonon spectral density^, a^F{ijj), namely: 

A = 2 /" "'^'"'"^ (25) 

where 
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a^F{uj) = p'^°\i2) ■ \U\ ■ -Im [D{uj)] (26) 

TT 

The real-axis form of the phonon propagator, D{uj), is calculated from its imaginary-axis values, 
Diiujy), by a Fade analytic continuation. 



III. RESULTS 

In choosing the parameters used to carry out the calculations, a number of criterion had to be 
satisfied. Firstly, we wished to operate outside the weak-coupling regime, so that vertex corrections 
would not be negligible. The phonon frequency needed to be large compared to conventional low- 
temperature superconductors, but not so large that it gave no realistic point of contact with those 
superconductors mentioned in the introduction. So we choose VL = equal to one-twelfth of the 
bandwidth. We had to ensure the electron-phonon coupling strength was not so strong that the 
ground state would contain bipolaronsS, making the perturbation expansion about a Fermi liquid 
state invalid. A maximum coupling oiU = —2t, hence a bare A < 0.5 ensured this. Finally, in order 
to calculate properties in an achievable time, while ensuring the imaginary frequency cut-offs were 
at energies larger than the band-width, the temperature of the calculations could not be too small, 
hence T > 2 x lO^^t. The last condition meant that results within the superconducting state were 
best carried out for as large U and fl as possible, so that Tc would be high. Hence, calculations 
near the band-edges, where the density of states was low, prove to be inaccurate, due to the very 
low transition temperatures there. 

Our first result is that dressing the 'bare' phonon propagator leads to considerable renormal- 
ization effects. To be specific, the value of A doubles from it's 'bare' value of A = 0.21, to A ~ 0.4 
after dressing the phonons, using parameters Q = t and U = — 1.5t, near half- filling. Moreover, at 
the increased interaction strength of U = —2t, A is enhanced by a factor of three from the value of 
A = 0.28 for undressed phonons to the dressed value of A ~ 0.9. Such an enhancement indicates 
that the perturbation expansion would be inaccurate at bare coupling strengths lower than might 
be naively expected. Figure |^ shows how a'^F{uj) is altered from its 'bare' value, a delta-function 
situated at uj = Q = t, when it is dressed. Note that there is both a shift to lower frequencies as 
well as a broadening of the spectrum. The shift to lower frequencies shows that a Holstein model 
with 'bare' phonon frequencies of near 10% of the bandwidth can be required to lead to dressed 
phonon frequencies at approximately 5% of the bandwidth. Hence the exact phonon self-energy 
used is an important factor when modeling electron-phonon systems, near the crossover between 
the weak-coupling and strong-coupling regimes. 

The preceding paragraph explains some of the dramatic differences between the critical tem- 
perature (Tc) values for the different approximations, shown in Figure |^. In particular, the Tc for 
dressed phonons is markedly higher than that for bare phonons, which is to be expected as A is 
also higher. Note that the transition temperatures fall rapidly with increasing filling, above a fill- 
ing of about n = 1.5, as the electronic density of states drops significantly in this region. Vertex 
corrections seriously reduce Tc near half- filling [n = 1). The two curves with dressed phonons 
in Figure |^ show a greater disparity than the two curves with 'bare' phonons in the same figure, 
which means that dressing the phonons, which increases the effective coupling, enhances the effect 
of vertex corrections. Although it is hard to distinguish the curves due to the low Tc near the band 
edge, above a filling of n = 1.75 the vertex corrections do lead to an enhancement of T^. This result 
is in agreement with previous worklllO. Note that all results show particle-hole symmetry, that is, 
they are symmetric about half-filling. The figures just show half of the band (n > 1), neglecting a 
mirror image below half-filling. 
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It is clear that vertex corrections do change Tc by a considerable amount, but for any experimen- 
tal measurement, where the microscopic parameters are not known, a Coulomb pseudopotential, 
IIq, can always be fitted to give the same Tc as vertex-corrections. Hence we continue with other 
results, to see where vertex corrections can not simply be mimicked by an appropriate /i^, which 
would cause the effects of vertex corrections to be unobservable. So we fit a fx^ to give the same 
value of Tc as vertex corrections, and go on to change the unobservable electron-phonon coupling 
strength, U, to give the same A as vertex corrections, when phonons are dressed. Hence the effects 
of vertex corrections can show up as discrepancies over a range of quantities compared to the values 
obtained with a fitted fi^ and A. 

Figure ^ is a direct comparison between the effects of vertex corrections and a Coulomb pseu- 
dopotential, /i^, on the value of the gap parameter. Our first method is to fix the bare electron- 
phonon coupling and to adjust the value of Uc and hence /i^ in a calculation without vertex 
corrections, until the same Tc is reached as found in the calculation with vertex corrections (shown 
in Figure 0(b)). Uc is then used unchanged, to calculate other properties such as the gap parameter. 
The Coulomb pseudopotential /i^ varies with electron filling at Tc as shown in the inset. Notice 
that an artificial value of fi^ < is required when n > 1.7, as vertex corrections enhance Tc in 
this region. Near half-filling, where fi^ is positive, and reduces Tc like vertex corrections, the gap 
parameter is reduced by a smaller amount. Hence a Coulomb repulsion leads to slightly higher gap 
ratio than vertex corrections. 

The second curve, with a lower value of /i^, is an alternate approach, where both the dressed 
value of A and Tc are fitted to the results with vertex corrections. A simple fit to Tc with a fixed 
bare coupling leads to a higher A with a Coulomb repulsion than with vertex corrections, because 
A is determined through the dressed phonon propagator. As a conventional analysis would fit A to 
the experimental data, as well as /i^, this second method follows the spirit of our paper by trying 
to fit conventional theory to the vertex-corrected results. In order to give the same values of A, 
the electron-phonon interaction energy had to be varied, and in fact reduced by 10% at half-filling. 
The result is a lower value of the gap ratio than with only fi^ fitted, but still a slightly larger value 
than with vertex corrections alone. 

Although the magnitude of the gap varies considerably, depending upon the approximation used. 
Figure |^ shows that the gap ratio, 2A/kTc, varies less markedly. The gap ratio is greater than 4 
in the case of dressed phonons without vertex corrections, which is typical of the strong coupling 
regime (A > 0.5). Note that when the phonons are 'bare', so the coupling is less strong, the inclusion 
of vertex corrections, while strongly reducing Tc and A individually, has little effect on the ratio, 
2A/fcTc. 

The isotope coefficient, a, has a value of 0.5 in the simplest, BCS, approximation, and in Migdal- 
Eliashberg theory with no Coulomb repulsion. The reason is that the phonon frequency provides 
the only cut-off for the coupling between different states, and phonon frequencies are proportional 
to M~°'^, where M is the ionic mass. Inclusion of a frequency-independent Coulomb repulsion, Uc, 
leads to a reduction in a, as does a finite bandwidth. The reduction in a, indicates that the increase 
in phonon frequency is less effective at increasing Tc than otherwise. Higher-frequency phonons re- 
duce the retardation in the electron-electron attraction, so the Coulomb repulsion between electrons 
is less shielded. The finite bandwidth means that the number of states coupled together includes a 
factor independent of phonon frequency, so Tc does not increase with Q as it might if there were an 
infinite number of electron states extending through all energies. 

It is known§'§ that Migdal-Eliashberg theory with a finite bandwidth and including a Coulomb 
repulsion, leads to the allowed rangjs of values for the isotope coefficient, < a < 0.5. Including 
a non-constant density of states^S can in principle lead to any positive value of a. The reason 
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being that Tc increases because extra electron states near the chemical potential are able to couple 
together when the phonon frequency increases. If the density of states increases significantly in 
the region of energy where new states are coupled together, the increase in Tc is much higher than 
would otherwise be expected, and a can be large, even greater than 0.5. The corollary is that if 
the density of states decreases away from the chemical potential, a also decreases, but never to less 
than zero, as T^. does not go down when the number of states coupled together goes up. 

Figure |10] shows that the inclusion of vertex corrections in the calculations with dressed phonons 
not only reduces a, but can in fact lead to negative values. Indeed, the strongest reduction in a by 
vertex corrections occurs near half-filling, and at strong coupling, where is comparatively large. 
By comparison, in all cases the Coulomb pseudopotential, which gives the same reduction in 
as vertex corrections, leads to a much smaller reduction in a. In Migdal-Eliashberg theory, a very 
small value of a requires a very low Tc. 

Hence, any observation of isotope effects which have a < 0, or a small a with moderate to 
high Tc, implies that either vertex corrections are involved or some other mechanism outside of 
Migdal-Eliashberg theory is important. Paramagnetic impurities0'0, proximity effect^!, anhar- 
monicity@'El, and an isotopic dependence of the electron density in the conduction bancoS, can 
also lead to a low or negative a without requiring a low T^. One or more of these effects may be 
important, in those materials with anomalously low values of but vertex corrections should 

also be considered. 

An important effect of dressing the phonons is that a small increase in the bare phonon frequency 
does not just result in a constant shift of o?F{ijj\ through a rescaling of the frequency variable. 
When calculating the isotope effect, it is common to assume that a mass substitution simply rescales 
the frequency, otherwise maintaining the same form of o?F[ili). However, the phonon self-energy is 
not independent of frequency, so the magnitude of o?F{uj) at its peak, which is inversely proportional 
to the imaginary part of the self-energy at that frequency, does not remain constant. In fact, the 
peak height is reduced by an increase in peak frequency, resulting in a slight reduction in A. Hence 
the increase in Tc with bare frequency is less than otherwise expected, reducing the isotope coefficient 
for dressed phonons. 

Interestingly, when the phonons are undressed, the Coulomb repulsion leads to a very small 
increase in a. This arises, because the term with a single polarization bubble (the first order term 
coming from dressed phonons) is present. The term acts to reduce a, but is less significant at the 
lower transition temperatures caused by the Coulomb pseudopotential. Meanwhile, the increase in 
yU* with temperature, which acts to reduce a is a much smaller effect. 



The free-energy difference, AT = Fg — Fn, is plotted as a function of filling, n, in Figure 11. In 



the simplest picture, true in the weak-coupling limit, one expects the magnitude of the free-energy 
difference to be approximately equal to p(/x)A^/2, representing a number of states, p(/i)A, each 
shifted by an average energy of order A/2. Although, with given by Z(0)p(°)(/i), the weak- 
coupling result predicts too high a condensation energyil, it does explain the qualitative changes 
between the different curves of Figure 0. 

In fact, the dimensionless quantity, '~fT^/{87iAF) [the Sommerfeld constant, 7 = 
27r^/c^p*^'')(/x)Z(0)/3] changes little for these curves. At half-filling, with dressed phonons, the value 
is reduced from the BCS constant result of 0.168 to a strong-coupling value of 0.137. The value 
with vertex corrections is 0.157, while with the Coulomb pseudopotential it is 0.149. Hence, as with 
the gap ratio, near half-filling, vertex corrections cause quantities to be closer to the weak-coupling 
values than does a Coulomb pseudopotential fitted for the same Tc. The values at a filling of n = 1.6 
are all closer to the weak-coupling limit, as expected when the density of states falls. The result 
is 0.163 for dressed phonons without vertex corrections, changing little to 0.165 with a Coulomb 
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pseudopotential and 0.158 with vertex corrections. 

It is worthwhile pointing out that the difference in thermodynamic potentials, which is usually 
calculated as an approximation to the free-energy difference, leads to very different results away from 
half-filling. The approximation is based on the assumption that the chemical potential changes little 
between normal and superconducting states, but this is not necessarily the case when there is a non- 
constant density of states. In fact, there is a particularly large shift from the normal-state chemical 
potential, /iTv, to that in the superconducting state, /i5, if /i^v lies near the van Hove singularity, 
where the non-interacting density of states is falling precipitously. This can be understood, by 
considering how states above and below the normal-state chemical potential, fi^, couple together 
and create an energy gap. When the density of states is falling rapidly with increasing energy, 
electronic states from a larger energy region above /iat couple to those in a small region below /i^v. 
The resulting energy gap is skewed up in energy, so the chemical potential in the superconducting 
state, Us, which sits in the middle of the gap, becomes greater than fiN. In the mirrored example 
below half-filling, near n = 0.4, f^s is also pushed away from half-filling, so we find < f^N there. 

The thermodynamic critical field. He, is effectively the square root of the free-energy difference, 
so shows qualitatively the same effects. is knowncZl to vary with temperature in a manner close to 
the behavior Hc{T) = Hc{0) [1 — (T/T^)^] which corresponds to the two-fluid model. The deviation 
function, plotted in Figure |I2], is the difference between the reduced critical field, Hc{T) / Hc{0) , and 
the quadratic fit, 1 - {T/T^f. 

The curve without vertex corrections at half-filling shows a small positive deviation, with a 
maximum of 0.02, typical of strong-coupling superconductors. Interestingly, when vertex corrections 
are included, the deviation function at half-filling becomes negative, showing a minimum of about 
—0.03, which is more typical of weak-coupling superconductors (BCS theory predicts a minimum of 
—0.037). This result fits in with those for other properties, demonstrating that vertex corrections 
reduce the effective coupling strength, near half-filling. A Coulomb pseudopotential, with the same 
power to reduce as vertex corrections, also decreases the deviation function, but to a lesser extent 
than vertex corrections do so. This is still true when the value of A is also fitted, by altering U as 
necessary. The calculations away from the band-center, at n = 1.6, lead to a negative deviation for 
all approximations. The reduced density of states at this filling leads to weak-couphng behavior. 

Other thermodynamic quantities, which can be derived from the free-energy data, are affected 
in similar ways. That is, vertex corrections reduce the effective coupling strength, to a greater 
extent than does a Coulomb repulsion giving the same T^. For example, vertex corrections reduce 
the specific heat jump, AC at T^., as does a Coulomb repulsion to a lesser extent. The following 
results are obtained by a numerical differentiation, so are not completely accurate in themselves, 
(perhaps only to 10%) but as much of the error is systematic, the trends are reliable. At half-filling, 
the dimensionless quantity, I\C j^T^ is reduced from the strong-coupling value of 2.44, to 1.63 with 
vertex corrections and 1.88 with a Coulomb pseudopotential (the BCS result is 1.43). Again, notice 
the typical result that near half-filling, vertex-corrections lead to a weaker-coupling result than does 
inclusion of a Coulomb pseudopotential. At an electron filling of n = 1.6, the results indicate less 
strong coupling, giving AC/7TC = 1.66 for dressed phonons, with the value reduced to 1.53 by 
vertex corrections, and to 1.44 by a Coulomb repulsion. 

IV. CONCLUSIONS 

We have completed a numerical investigation of the effects of vertex corrections, dressing phonons 
and a non-constant density of states on the physical properties of strong-couphng superconductors. 
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We solved the Holstein model, using four distinct perturbation theories, within conserving approx- 
imations. 

The necessity of incorporating a realistic phonon self-energy is of considerable importance to 
those working with model Hamiltonians. The use of dressed phonons in the Holstein model, leads 
to a large renormalization of the parameters - in particular, the value of A can be enhanced by a 
factor of 3, when its 'bare' value of 0.28 would suggest the system is in the weak-coupling regime. 
Such an enhancement of A reveals itself in increased T^, a, A, He, and a gap ratio (2A/fcTc) greater 
than 4. The real-frequency data shows that the Einstein spectrum (a delta-function at Q) is both 
broadened and peaked at a lower frequency, when the bare, Einstein phonons are dressed. 

The non-constant density of states affects Migdal-Eliashberg results in both expected and un- 
expected ways. Firstly, all quantities which depend on the density of states as a parameter within 
Migdal-Eliashberg theory change in the expected manner as the electron band-filling changes. Note 
that any sharp features in the normal electronic density of states have their effects reduced by the 
'averaging' over a large phonon frequency range. Hence the strong fall in Tc and A is expected at 
both small and large fillings {n < 0.45 and n > 1.55) due to the rapid fall in the density of states in 
our model. More subtle, is the result that the chemical potential shifts by a considerable amount 
between the normal and superconducting states, if it lies near a van Hove singularity. To observe 
such an effect, the superconductor would have to be coupled to one with a more constant density 
of states. 

We find, in agreement with previous work0'001ii, that vertex corrections lead to results that 
correspond to a reduced effective strength of the electron-phonon coupling near half-filling, but 
an increased coupling strength near the band edges. These effects are exemplified by reductions in 
critical temperature, Tc, superconducting gap. A, isotope coefficient, a, and thermodynamic critical 
field. He, near half-filling. As nearly all of these effects can be modeled by an appropriate Coulomb 
pseudopotential, /ij., it makes it extremely difficult for any single experiment to reveal that vertex 
corrections have played a significant role. However, we do find some trends worth pointing out. 

Firstly, if there is difficulty in fitting both Tc and A with a given a'^F{uj) and /i^, then this is 
an indication that vertex corrections may contribute, since they affect the ratio 2A/kTc for fixed 
Tc, reducing it near half-filling. Secondly, if the experimentally measured deviation function for the 
thermodynamic critical field lies below the predicted value [with a given a'^F{uj) and /i^] vertex 
corrections could be important. Thirdly, a theoretical prediction, ignoring vertex corrections, will 
overestimate the specific heat jump at Tc. Finally, a more striking result, is in the isotope coefficient, 
a, which vertex corrections reduce much more markedly than does /i^. Indeed, vertex corrections 
can lead to a < 0, which /i^ alone can never dolll. Hence materials that have moderate to large TcS, 
but small isotope coefficients can still be electron-phonon mediated superconductors with vertex 
corrections included. As the isotope coefficient is the single experimental quantity affected the 
most by vertex corrections, it is important to consider materials which have a unexplained by 
Migdal-Eliashberg theory. 

Anomalously low, and even negative isotope coefficients have been observed in materials, such as 
Ru0, a-uranium0, Pdni, and La2-:,Sr^Cu04ii^, where vertex corrections are probably not im- 
portant, and other mechanisms, such as anharmonicity, conduction electron density variations and 
paramagnetic impurities play a role. However, a system such as RbaCeol^^, where the phonon fre- 
quency is a sizable fraction of the electron bandwidth, is much more likely to have vertex corrections 
affect the value of a. 

Still, the best way to see the effects of vertex corrections is to directly view their contribution in 
the multiphonon region of a tunnel junction, or in the high-energy region of the optical conductivity. 
So, in order to unequivocally demonstrate the presence of strong vertex corrections in a material. 
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more accurate dynamical measurements at energies beyond the highest phonon frequencies need to 
be carried out. 
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APPENDIX A: CONSERVING APPROXIMATION FORMULA 



In this appendix, we give the specific formulae for the free-energy functional, self-energy and 
irreducible vertex functions for each of the four conserving approximations. Figures ^-|^ are the 
representations of these equations as Feynman diagrams. Hereafter, we employ the shortened 
notation G„ = G{iuin), G*^ = G{—iuJn) and similarly for F„, D,^ and vr^^-'. Note that the difference of 
two fermionic frequencies leads to a bosonic frequency, i.e. Dm-n = D{iuJm — i^^n) = D{iujy) where 
V = m — n. 

The calculations with no vertex corrections and a 'bare' phonon propagator have the free-energy 
functional, 



bare 



(Al) 



VTrlr^G T^G ]D. 



(0) 



where vr^^-' is the electron polarizability, defined in Eq. (H). The inclusion of the matrices, ensures 
that each pair of off-diagonal Green's functions, F*Fm, corresponding to Cooper pair creation then 
annihilation, enters the product with a minus sign. 

Functional differentiation with respect to the diagonal electron Green's function, G„, and off- 
diagonal Green's function, F* leads respectively to the diagonal term in the self-energy, E(ic<j„): 



D 



(0) 



(A2) 



and the off-diagonal term, 

ibare / 



UT E FrnD^^l^ - U^T E F^Til^ 



(0) 



D 



(0) 



(A3) 



We only require the superconducting vertex part, which is given by the derivative: 

TF^';- = 6<P{zu^)/6F^, (A4) 

taken in the limit Fm ^ 0. Hence with bare phonons and no vertex corrections, the vertex function 
is that shown in Figure ^(a): 



■^bare 



D 



(0) 



(A5) 
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There is only one extra term which comes from the inclusion of vertex corrections in the free- 
energy functional. It is the diagram in Figure ^(b) with crossed phonon lines, and is equal to: 



4 



n,m,l 



T G T G T G, T G , , 



r,{o) ^(0) 



(A6) 



The total free-energy functional for bare phonons with vertex corrections included is + (^'"^. 

The extra term in the self-energy is found by differentiation of the above term. For the diagonal 
and off-diagonal parts, respectively, this leads to: 



E-(^a;„) = [/2T2^{a 



ml 



1^(0) ^(0) 
J -'-^n—m-'-^ 



(A7) 



\ TT2rp2 \ ^ { rp \ rp* rp H H 

m.l 

-F/GmC^-m+i} (A8) 

The extra term leads to the new self-energy, S = S''"'"^ + S^^ and = 0^°^^*^ + 0"'^. 

When the vertex correction term is added to the self-energy, three new terms appear in the 
irreducible vertex function, coming from each of three Green's functions that can be differentiated. 
The extra diagrams in Figure ||(b) contribute a total of: 

^m,n = U'^T'Y I \^—GiGn-m+i — G*^_^^iG*i Dnl^D\}^ — G iG*i^^j^_iD[^}_iD\J_^ , (A9) 

leading to Tm,n = ^m^n + ^mn irrcduciblc vertex function for 'bare' phonons, with vertex 

corrections. 

The calculations with no vertex corrections, but a dressed phonon propagator have the free- 
energy functional, 

^dressed = ^— ^ Tr gjTr £ J/^(°) (o; = 0) (AlO) 



2 
UT 



n,m 

2 



VTrfG G ]Dr, 



Functional differentiation with respect to the diagonal electron Green's function, G„, and off- 
diagonal Green's function, leads respectively to the diagonal term in the electron self-energy, 

^dressed/- V 



^dressed. 



Un + UTY,G.^D,,^^ (All) 



and the off-diagonal term, (j)'^'^'^'^^^'^{iujn)'- 

^dressed^^^^^ = -UT FmD„.-m. (A12) 

m 

Similarly, differentiation with respect to the dressed phonon propagator leads to the phonon self- 
energy. 
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n(^) {luj,) = -2UTY^ [GmGm+u - J ( A13) 

m 

where the factor of 2 indicates a sum over electron spins. 

The superconducting vertex function still retains its simple form, differentiation of Eq. ( [A12| ) 
giving: 

^dressed ^ _JJD^_^. (Al4) 

Analogously to the case of 'bare' phonons, there is only one extra term which comes from vertex 
corrections in the free-energy functional. It is equal to: 

= ^ E Tr [r3 a„ 13 13 Is a„-™+J D^_:^D,_^, (A15) 

n,m,l 

where now the total free-energy functional for dressed phonons with vertex corrections included is 

(^dressed _j_ ^fc2 

Differentiation of the above contribution now leads to an extra term in the phonon self-energy, 
n'''^, as well as the extra electronic self-energy terms. 11^'^ is the term shown with crossed phonon 
lines in Figure |5(b): 

U-i^u^) = -U-T- E Tr {r, g^^^ Z, Z, Z, g^J D,_^. (A16) 

m,,l 

Hence the full phonon self-energy is now Ii{iuj^) = Il^^\iuj,y) + IF^{iuj^). 

The vertex functions and self-energies for electrons gain terms equivalent to those in Eqs. ( [A7| - 
|A9| ) only with bare phonon propagators, -D^°\ replaced by dressed ones, D^. 

In the normal state, the expressions for the free-energy functionals and self-energies are simpli- 
fied, by setting all off-diagonal contributions to zero, = and (piiujn) = 0. Note that the above 
formulae for irreducible vertex functions are only calculated in the normal state. 
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FIGURES 
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FIG. 1. The single-spin density of states for non-interacting electrons on a 3D cubic tight-binding 
lattice. Note the nature of the van Hove singularities, which lead to an abrupt fall in the density of states. 
They occur at electron fillings of n ~ 0.45 and n ~ 1.55. The approximate form used in the calculations is 
also plotted, and is indistinguishable by eye from the exact curve. 
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a) 



b) 



c) 



d) 




FIG. 2. Different forms for the electron self-energy. All diagrams include the same initial Hartree 
term, followed by the specific Fock contribution and any other required diagrams. All electronic Green's 
functions are dressed, but phonon lines can be undressed (thin) or dressed (thick) as in Figure |,b). For 
bare phonons, (a) is with no vertex corrections, while (b) includes the vertex-corrected diagram. When 
the phonons are dressed, (c) is without vertex corrections, while (d) includes them. A factor of 2 appears 
in front of each electron loop term, to indicate a sum over spins. 
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^AAAf =^AAA/ ^AAAJ(^^)^AAA/ 



b) 



<5>-0-<D 



FIG. 3. Dyson's equation for the phonon propagators. Thin hnes indicate bare Green's functions while 
thick hnes indicate dressed ones, (a) is the phonon self-energy without vertex corrections, (b) is with vertex 
corrections. The factors of 2 come from a sum over electron spins. 
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a) 



b) 



c) 



d) 




FIG. 4. Free-energy functional, $ for the conserving approximations. Equations (a) and (b) are with 
bare phonon propagators, while (c) and (d) have dressed phonons. The first term on the right in each 
equation is the Hartree term, which does not change and which is not included in (Its phonon line is 
always undressed, to avoid double counting). The extra, final diagram in (b) and (d) is the vertex-correction 
term. 
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FIG. 5. The irreducible vertex-function diagrams, for superconductivity, (a) is with no vertex correc- 
tions, while (b) includes the three vertex-corrected diagrams. When the phonons are dressed, (c) is without 
vertex corrections, while (d) includes them. 
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0.60 




Frequency, CO / ^ 

FIG. 6. a^F{u}), calculated from the dressed plionon propagator on the real frequency axis for n = 1 
and Q = t. Increased coupling leads to a greater downwards shift in frequency from the initial delta-function 
at Lo = t. 
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FIG. 7. Transition temperature as a function of filling. The dashed lines indicate vertex corrections are 
included, while solid lines are without them. The lines with circles have dressed phonons. Both diagrams 
show an enhanced Tc by dressing the phonons, while vertex corrections reduce Tc near half-filling. In all 
cases, Tc falls rapidly above a filling of n = 1.55, where the density of states drops rapidly. 
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FIG. 8. A comparison between the effects of vertex corrections and a Coulomb pseudopotential, fi^ 
on the superconducting gap. The calculations arc with dressed phonons, with U = —2t and = t. The 
dashed curve is with vertex corrections, while the dotted curves include a /z^, whose value changes with 
filling as shown in the inset, to ensure the two corresponding curves are exactly the same. The dotted 
curve with triangles indicates a fit to the same A by adjusting U as well as jiQ. 
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FIG. 9. Gap ratio, 2A/kTc, as a function of filling, n. Dressed phonons (with circles) exhibit strong 
coupling behavior, by the increased gap ratio. Vertex corrections (dashed lines) reduce the effective coupling 
strength in the center of the band, hence the gap ratio is lower in this region. The dotted lines indicate 
how a Coulomb pseudopotential alters the gap ratio, (when Tc is matched to Tc with vertex corrections, 
and where triangles indicate that A is also matched, by adjusting U). With dressed phonons, the Coulomb 
repulsion clearly has less of an effect than do vertex corrections, but this is not the case with bare phonons. 
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FIG. 10. The isotope coefficient, a, plotted against electron filling, n, for (a) U = —1.5t and (b) 
U = —2t. Vertex corrections, indicated by dashed lines, reduce a. As does dressing the phonons, seen by 
the lines with circles. Figure (b) exhibits the unusual feature of cc < for dressed phonons, with vertex 
corrections included. The dotted curve in (b) is the result with the Coulomb repulsion, /x^, shown in the 
inset in the previous figure. 
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FIG. 11. Gondensation energy, Fs — F^. Circles indicate the phonons are dressed, which leads to a 
greater condensation energy. The dotted lines indicate dressed phonons with a Coulomb pseudopotential, to 
mimic the vertex-corrected Tc, the triangles indicating that U is also adjusted to mimic the A obtained with 
vertex corrections. The dashed, vertex-corrected curve shows smaller condensation energy for n < 1.55, 
though the difference can not be seen with bare phonons. 
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FIG. 12. Critical-field deviation function, Hc{T) / Hc{0) — [1 — (T/Tc)^]. Solid lines are without, dashed 
lines are with vertex corrections. The dotted lines with and without triangles include the Coulomb pseu- 
dopotentials, /Li^ of Figure ^ inset. All results are for dressed phonons. (a) is at half-filling, n = 1 while 
(b) is at a filling of n = 1.6. Note the shift in scale for the second graph. 



28 



